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From a combination of probes including the cosmic microwave background (WMAP7+SPT), 
Hubble constant (HST), baryon acoustic oscillations (SDSS+2dFGRS), and supernova distances 
(Union2), we have explored the extent to which the constraints on the effective number of neutrinos 
and sum of neutrino masses are affected by our ignorance of other cosmological parameters, in- 
cluding the curvature of the universe, running of the spectral index, primordial helium abundance, 
evolving late-time dark energy, and early dark energy. In a combined analysis of the effective num- 
ber of neutrinos and sum of neutrino masses, we find mild (2.2cr) preference for additional light 
degrees of freedom. However, the effective number of neutrinos is consistent with the canonical 
expectation of 3 massive neutrinos and no extra relativistic species to within la when allowing for 
evolving dark energy and relaxing the strong inflation prior on the curvature and running. The 
agreement improves with the possibility of an early dark energy component, itself constrained to 
be less than 5% of the critical density (95% CL) in our expanded parameter space. In extensions 
of the standard cosmological model, the derived amplitude of linear matter fluctuations a$ is found 
to closely agree with low-redshift cluster abundance measurements. The sum of neutrino masses is 
robust to assumptions of the effective number of neutrinos, late-time dark energy, curvature, and 
running at the level of 1.2 eV (95% CL). The upper bound degrades to 2.0 eV (95% CL) when 
further including the early dark energy density and primordial helium abundance as additional free 
parameters. Even in extended cosmological parameter spaces, Planck alone could determine the 
possible existence of extra relativistic species at 4<j confidence and constrain the sum of neutrino 
masses to 0.2 eV (68% CL). 



I. INTRODUCTION 

Observations of the cosmic microwave background 
(CMB) [3-11], large-scale structure and type la 

supcrnovac (SNe) have established a flat ACDM 
model, with nearly scale-invariant, adiabatic, Gaus- 
sian primordial fluctuations as providing a consistent 
description of the global properties of our universe. 
At the same time, we do not yet understand the mi- 
croscopic identities of the dark energy (A), cold dark 
matter (CDM), and inflaton (primordial fluctuations) 
that enter our standard cosmological model. 

The neutrino sector is another area that the stan- 
dard model is yet unable to fully describe, with open 
questions related to the effective number of neutrinos 
N e s and their masses m v . A joint analysis of CMB 
data from WMAP7 with baryon acoustic oscillation 
(BAO) distances from SDSS+2dF and Hubble con- 
stant from HST reveals a weak preference for extra 
relativistic species (N cS = 4.34 ±0.87) 0. When fur- 
ther combined with small-scale CMB data from ACT 
or SPT, this preference mildly increases to the 2a level 
(N e3 = 4.56 ± 0.75 with addition of ACT [l| and 
N eS = 3.86 ± 0.42 with addition of SPT [ll|; further 
see [T^ - [27{ ) . A primary objective of this manuscript 
is to clarify how robust these recent indications of ad- 
ditional light degrees of freedom are to assumptions 
of the underlying cosmology, in particular to alterna- 
tive models of the dark energy, curvature of the uni- 
verse, running of the spectral index, primordial helium 
abundance, and to the sum of neutrino masses, which 
we know is nonzero from neutrino oscillation experi- 



Parameter 


Symbol 


Prior 


Baryon density 


Q b h 2 


0.005 


-> 0.1 


Cold dark matter density 


Q. c h 2 


0.01 - 


> 0.99 


Angular size of sound horizon 


e s 


0.5 - 


■> 10 


Optical depth to reionization 


T 


0.01 - 


■> 0.8 


Scalar spectral index 


n s 


0.5 - 


> 1.5 


Amplitude of scalar spectrum 


ln(10 10 ^) 


2.7 - 


->■ 4 


Effective number of neutrinos 


iVeff 


1.047 


->■ 10 


Sum of neutrino masses 


Em, [eV] 


- 


> 5 


Constant dark energy EOS 


w 


-3 





Running of the spectral index 


dn s 
dink 


-0.2 


-)• 0.2 


Curvature of the universe 


n k 


-0.4 - 


-+ 0.4 


Primordial helium abundance 


Y P 


- 


> 1 


Present dark energy EOS 


w 


-3 - 


-+ 


Derivative of dark energy EOS 


W a 


-10 


10 


Early dark energy density 


fie 


-> 


0.2 



TABLE I. We impose uniform priors on the above cosmo- 
logical parameters. In addition, we always consider the 
Poisson point source power D^oo? the clustered power 
D3060! and the SZ power -Dfooo as nuisance parameters 
constrained by the CMB data [TJ|. Moreover, we always 
derive as, the amplitude of linear matter fluctuations on 
scales of 8 Mpc/h at 2 = 0. Beyond a constant dark energy 
equation of state (EOS), we also consider a time- varying 
expansion w(a) = wo + (1 — a)w a , and an early dark en- 
ergy model described in Sec. IIII Dl In this table, the first 
6 parameters are defined as "vanilla" parameters. 



ments |28l-l31|. 

Constraining N e g with cosmology is mainly 
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TABLE II. Constraints on Cosmological Parameters using SPT+WMAP+Hq+BAO. 







ACDM 


ACDM 


ACDM 


ACDM 


ACDM 


ACDM 








+ AW 


+ E m f 




+N cft +Y p ■ 


fA r off±E m ^± y p 


± i iiiicLi y 


1 nflOi h' 2 


2 237 + n 038 

— .—,■>) _1_ U.UJO 


2 261 + 042 


2 238 + 039 


2 272 + 043 


2 271 + 043 


2 273 + 044 




i noo /j 2 


ii 99 _l n oe 




ii ii -u n 9Q 


i q nc; _|_ n Q4 


12.5 ± 1.1 


13.0 ± 1.2 




10 4 6» s 


104.12 ±0.15 


103.95 ± 0.17 


104.15 ± 0.15 


103.95 ±0.18 


104.07 ±0.29 


103.96 ±0.29 




r 


0.086 ± 0.014 


0.086 ±0.014 


0.088 ±0.014 


0.090 ±0.015 


0.088 ±0.014 


0.090 ±0.015 






0.9648 ± 0.0092 


0.981 ±0.013 


0.9661 ± 0.0096 


0.987 ±0.013 


0.983 ±0.013 


0.987 ±0.013 




^(lO 10 ^) 


3.195 ± 0.034 


3.186 ±0.035 


3.184 ±0.035 


3.170 ±0.037 


3.180 ±0.036 


3.169 ±0.037 


Extended 


AW 




3.87 ±0.42 




4.00 ± 0.43 


3.70 ±0.54 


3.99 ±0.59 




"Emu [eV] 






< 0.45 


< 0.67 




< 0.73 




y P 










0.277 ±0.037 


0.261 ±0.039 


Derived 


°"8 


0.811 ± 0.018 


0.862 ±0.033 


0.758 ±0.042 


0.798 ±0.053 


0.860 ± 0.034 


0.796 ±0.055 



Mean of the posterior distribution of cosmological parameters along with the symmetric 68% confidence interval about the 
mean. We report the 95% upper limit on the sum of neutrino masses E vn v . The primordial helium mass fraction Y v is enforced 
consistent with standard BBN unless we allow it to vary as a free parameter. 



achieved through tight CMB measurements of the red- 
shift to matter-radiation equality z eq , the baryon den- 
sity Q.bh 2 , the angular size of the sound horizon 6 S , and 
the angular scale of photon diffusion 9d [HJ • Keeping 
z eq and fifth 2 fixed as N e g increases can be achieved 
by increasing the dark matter density Q c h 2 (assum- 
ing massless neutrinos), which manifests in a large 
correlation with N e g (shown in Fig. [5]). Meanwhile, 
an increase in N e g and Yo both yield an enhanced 
Silk damping effect [H, [H-IH, and by fixing 9 S it 
can be shown that 6 d oc (1 + /„) ' 22 / y/l + Y p [j3], 
where /„ = p v / P~t is proportional to N c g. As a conse- 
quence, the suppression of the CMB damping tail can 
be picked out as a signature of extra relativistic species 
when Y p is known, while the constraints on N e g are 
relaxed when allowing for Y p as a free parameter. 

An increase in N e g further shifts the acoustic peak 
locations [34| . but this has been shown to be a small 
effect [l4j|. Instead, the constraint on N c g can be im- 
proved by the inclusion of low-redshift distances and 
a prior on the Hubble constant, H , as these are use- 
ful in constraining fl c h 2 and by extension N c g. How- 
ever, when allowing for evolving dark energy, the abil- 
ity to improve constraints on AT c ff from observations 
of the expansion history becomes diminished, as il- 
lustrated by the error ellipses for {N c g , Q c h 2 , w} in 
Fig. [2j Therefore, the inclusion of SN data becomes 
critical to a precise determination of the effective num- 
ber of neutrinos. 

The dark energy equation of state (EOS) is 
moreover anti-correlated with the sum of neutrino 
masses BUS! HI- In the CMB temperature power 
spectrum, the sum of neutrino masses shifts the first 
peak position to lower multipoles by changing the frac- 
tion of matter to radiation at decoupling, which can 
be compensated by a reduction in the Hubble con- 
stant (similar to the case for positive universal curva- 
ture) [1, |H, H3] ■ B AO distances and an H prior can 
therefore be used to reduce correlations between the 



sum of neutrino masses and the dark energy EOS, but 
also with the curvature density. 

The strongest limits on the sum of neutrino masses 
from the CMB combined with probes of the expansion 
history and matter power spectrum place it at sub-eV 
level [3-11, [37l - l53| . We take the conservative approach 
in only combining CMB data with low-redshift mea- 
surements of the expansion history. While SN ob- 
servations play an important role in constraining the 
dark energy EOS and thereby reduce the correlation 
between "^2 vn>u and w, these observations are not pow- 
erful in constraining the curvature of the universe and 
therefore less helpful in reducing the correlation be- 
tween "^2 m v and ■ 

Beyond the vanilla parameters and the three addi- 
tional parameters {-/V e ff , "^2 "V, w}, we relax the com- 
monly employed strong inflation prior on the univer- 
sal curvature Qk and running of the spectral index 
dn s /dhik. Given that most popular models of infla- 
tion predict \dn s /dlnk\ < 1CT 3 [H HU and |fi fc | < 
1CT 4 (e.g. [13, [M H3), at the level of precision of 
present CMB data it is generally justified to fix these 
two parameters to their fiducial values of zero. How- 
ever, given the mild preference for N e g > 3 0, [HI HH , 
we allow for the possible existence of inflationary mod- 
els with large curvature or running. In particular, 
~ I0~ 2 may be generated in models of open infla- 
tion in the context of string cosmology [54] . [58j , while a 
large negative running may be produced by multiple 
fields, temporary breakdown of slow-roll, or several 
distinct stages of inflation [54| l59M6l| . 

Among alternatives to the cosmological constant 
with w = — 1, the most popular are scalar field models 
with potentials tailored to give rise to late-time ac- 
celeration and current equation of state for the dark 
energy, w, close to -1 [62H69| . Like a cosmological con- 
stant, these models are fine-tuned to have dark energy 
dominate today. However, the requirement w > — I 
currently, does not imply that dark energy was sub- 
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TABLE III. Constraints on Cosmological Parameters using SPT±WMAP±# ±BAO. 







iuCDM 


ACDM 


«;CDM 


ACDM 


ibCDM 


wCDM 












±A r off±E m ^ 

± ± ^fc 
In k 


+N cff +J2m l , 
±7nrrr ± ^fc 

d In fc 


+N cft +Y,m„+Y p 

± JTTL ± 
d\n k 


JT 1 lllltll V 


lOOf^/i 2 


2.219 ±0.042 


2.272 ±0.043 


o 994 4. a nisi 


9 944 4- (1 (1^4 


9 1 09 4. n nfis 


9 1 f?a 4. n 07Q 

Z.1UO ZC U.Uit? 




ioon c /i 2 


11.44 ±0.45 


13.05 ±0.94 


1 2 96 + 96 


13.2 ± 1.0 


12.4 ± 1.2 


13.1 ± 1.7 




10 4 6» s 


104.09 ±0.16 


103.95 ±0.18 


1 (13 97 + 1 8 


1 nr> Q7 4- ri 1 9 


1 (14 (Ifi 4- (1 21 


1 01 fil + (1 44 




T 


0.083 ± 0.014 


0.090 ±0.015 


086 ±0015 


090 ± 015 


088 ± 015 


091 ± 016 




n a 


0.958 ± 0.011 


0.987 ±0.013 


n ors 4- n (122 


n 07a 4- n (11 B 

U.l/ 1 U _1_ U.U1U 


n qc;c; 4. n (125 


n 040 4. ri (127 




ln(10 10 A s ) 3.216 ±0.042 3.170 ±0.037 


3.211 ± 0.052 


3.179 ±0.045 


3.210 ±0.052 


3.200 ±0.057 


Extended 


w 


-1.10±0.11 




-1.31 ±0.30 




-1.46 ± 0.39 


-1.35 ±0.41 




N ee 




4.00 ±0.43 


3.59 ± 0.57 


3.74 ±0.58 


3.10 ±0.74 


3.38 ± 0.86 




Emu [eV] 




< 0.67 


< 1.2 


< 1.2 


< 1.2 


< 1.4 




dn s 
din fc 








-0.011 ±0.019 


-0.018 ± 0.019 


-0.033 ±0.031 




ioon fc 








0.75 ±0.93 


0.13 ±0.99 


0.76 ± 1.5 
















0.196 ± 0.084 


Derived 


CT8 


0.848 ± 0.049 


0.798 ±0.053 


0.775 ± 0.063 


0.768 ±0.070 


0.803 ± 0.085 


0.779 ± 0.091 



Same as Table ITlI but with the addition of {w, dn a /d In k, f2fc}. Due to the large correlation between n B and dn s /dlnk at our 
pivot scale ko = 0.002/Mpc, we quote values for n s at a less correlated scale ko = 0.015/Mpc. For the "uiCDM+A r e £f + 5~J m " ± 
dn s /dlnk + fij." case where N e g is closest to to the boundary at 3, we also considered a run where we impose a hard prior 
of N e ff > 3. Here, we find N c g = 3.65 3qq' 3 oo> wnere * ne * wo se * s °f upper and lower boundaries denote 68% and 95% CLs, 
respectively. The changes to the sum of neutrino masses and other parameters that weakly correlate with N c g are small (< 10%). 
While all within la, the largest changes are seen in 100f2 c /i 2 = 13.17±0.97 (compared to lOOf^/i 2 = 12.4±1.2), w = -1.25±0.30 
(compared to w = -1.46±0.39), n s = 0.971±0.019 (compared to n B = 0.955±0.025), dn a /d\ak = -0.0088±0.0168 (compared 
to dn B /dlnk = -0.018 ± 0.019), and 100f2 fe = 0.2 ± 1.1 (compared to 100C fe = 0.13 ± 0.99). This particular configuration of 
parameter space and datasets shows the largest extent to which parameters may change with an N e g > 3 prior as compared to 
our other runs. The changes to the parameters are more modest when including SNe because of the preference for larger values 
of N e g, as seen in Table Hvl 



dominant at earlier times, specifically redshift z > 2, 
where we have no direct constraints. Given the degen- 
eracy between dark energy and the sum of neutrino 
masses, we further consider a model that describes 
dark energy as non-negligible in the early universe in 
Sec-HTD] 

We describe our analysis method in Section 2. In 
Section 3, we provide constraints on a ACDM model 
with three massive neutrinos and additional light de- 
grees of freedom, then follow up with successive addi- 
tions of a constant dark energy equation of state, uni- 
versal curvature, running of the spectral index, and 
primordial helium abundance (all parameters defined 
in Table QJ . We also explore the constraints for a time- 
varying dark energy equation of state, including an 
early dark energy model. Lastly, we compare the con- 
straints from present data to that expected for Planck. 
Section 4 concludes with a discussion of our findings. 

II. METHODOLOGY 

We employed a modified version of CosmoMC [7(1 
|7"H in performing Markov Chain Monte Carlo 
(MCMC) analyses of extended parameter spaces with 
CMB data from WMAP7 Q and SPT [ll[, BAO dis- 
tance measurements from SDSS±2dFGRS [z3], the 
Hubble constant from HST [l2[, and SN distances 
from the SCP Union2 compilation [73[. All param- 



eters are defined in Table HI In determining the con- 
vergence of our chains, we used the Gclman and Rubin 
R statistic [74[, where R is defined as the variance of 
chain means divided by the mean of chain variances. 
To stop the runs, we generally required the conserva- 
tive limit (R — 1) < 10~ 2 , and checked that further 
exploration of the tails does not change our results. 

The CMB temperature and E-mode polarization 
power spectra were obtained from a modified version 
of the Boltzmann code CAMB [zl [7^. We approx- 
imated the effect of a dark energy component with 
a time-varying EOS by incorporating the PPF mod- 
ule by Fang, Hu, & Lewis (2008) [77] into CosmoMC. 
Given that the small scale CMB measurements of SPT 
come with much smaller error bars than ACT [ToL [Tl1 | , 
the further inclusion of the ACT dataset would not 
lead to significant improvements in our constraints, 
as we explicitly checked. 

When allowing for nonzero neutrino rest mass, 
we distribute the sum of neutrino masses m v = 
94 eV Q„h 2 ) equally among 3 active neutrinos. We 
treat additional contributions to N e g as massless, such 
that N e g = (3+iV m i), where iV m i denotes the massless 
degrees of freedom. Since we impose 1.047 < N e g < 
10, the number of rclativistic species is always posi- 
tive at early times. At late times, our prior on N e g 
implies that the number of relativistic species can be 
negative (—1.953 < N m \ < 7). However, the total ra- 
diation energy density (cx 1 + 0. 227-/V m i at late times) 
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TABLE IV. Constraints on Cosmological Parameters using SPT+WMAP+Ho+BAO+SNe. 







toCDM 


u?CDM 


wCDM 

1 n k 


ioCDM 
+N cft +Y,m lJ +Y p 
±77T^rr + ^fc 

d In A- 


L I llilell V 


1 DflOi./) 2 


o ooq _|_ n 041 


9 9K7 _l_ n rMQ 


o oofi -4- n n^iQ 


o 1 71 -L n 080 




ioon c ft 2 


ii sfi + n 41 


1 3 1 4 + 94 


13.1 ± 1.0 


14.0 ± 1.4 




1O 4 S 


104.10 ± 0.16 


103.95 ± 0.17 


103.99 ±0.18 


103.66 ± 0.36 






082 + 014 


088 + 015 


n naq + n 016 


090 + 015 




' s 


n qfin + 01 

U.C/UU _1_ U.U1U 


n n«i inmi; 

U. .701 _l_ U.uU 


n qvn -J- n 01 9 


+ n 026 




ln(10 10 A s ) 


3.209 ±0.039 


3.185 ± 0.041 


3.194 ±0.049 


3.197 + 0.053 


Extended 


10 


-1.049 ±0.072 


-1.09 ±0.11 


-1.10 ±0.11 


-1.13 ±0.12 




N cff 




3.88 ± 0.44 


3.58 ±0.60 


3.78 ± 0.61 




[eV] 




< 0.92 


< 1.2 


< 1.7 




dn s 
d In k 






-0.013 ±0.019 


-0.035 ±0.030 




won k 






0.64 ±0.95 


1.2 ± 1.1 












0.176 ±0.079 


Derived 


CT8 


0.830 ± 0.038 


0.790 ± 0.060 


0.774 ±0.072 


0.751 ± 0.081 



Same as Table IIIII but with the addition of supernova distance measurements from the Union2 compilation. For the case 
"uiCDM+A'cff +52 m^+dris/dlnfc+flfc," we also considered a run where we impose a hard prior of N e g > 3. Here, we find 
N e 3 = 3.74 3 qq' 3 oo> where the two sets of upper and lower boundaries denote 68% and 95% CLs, respectively. The largest 
changes this prior induces in other parameters are in n s = 0.974 ± 0.017 (compared to n s = 0.970 ± 0.019), dn s /d\nk = 
-0.010 ±0.017 (compared to dn a /d\nk = -0.013 ± 0.019), and 100C fe = 0.52 ± 0.92 (compared to 100C fe = 0.64 ±0.95). All 
other parameters are modestly affected by our choice of prior (< 10%). 



is always positive. We chose this particular prior on 
N e ff in order for the data itself to rule out a given 
part of parameter space. In Figs [2H1 we find that the 
marginalized contours on N c ff close before the lower 
end of our prior, such that the data itself is constrain- 
ing the radiation content from below. For complete- 
ness, we also considered several conventional runs with 
the prior N g > 3, such that N m \ > 0, and we find 
no qualitative changes in our results. For complete 
details, see the captions of Tables [TTTT IIV[ IVII 

As part of our analysis of extended parameter 
spaces, we consider cases with the primordial frac- 
tion of baryonic mass in helium Y p as an unknown 
parameter to be determined by the data. However, 
when we do not allow Y p to vary freely, it is deter- 
mined in a BBN-consistent manner within CAMB via 
the PArthENoPE code [zl], which enforces 

Y p = 0.2485 + 0.0016 [(273.9ft b /i 2 - 6) + 100(5- 1)] . 

(1) 

Here S = y/l + (7/43) AN„ encapsulates deviations 
from standard BBN [79t48ll ] . and we let AN V = (N eS - 
3.046) in agreement with the SPT analysis. Aside 
from the derived limits on Y pi we explicitly checked 
that our results do not significantly change (< 10%) 
when passing AN V = to PArthENoPE instead. 

Furthermore, in our analysis we cither consider "en- 
forcing the strong inflation prior" on the curvature 
and running, by which {Q^ = 0, dn s /d In k = 0}, 
or "relaxing the strong inflation prior" such that 
{fife, dn s /d\nk} are allowed to vary as free parame- 
ters to be constrained by the data. We define the run- 
ning of the spectral index via the dimensionless power 



spectrum of primordial curvature perturbations: 

/ h \ 7ls ~ ^ 2 l J *(k/ko)dn 3 /d In k 

A%(k) = A% ,{ko) (-J , (2) 

where the pivot scale ko = 0.002/Mpc. Due to the 
large correlation between n s and dn s /d\wk at this 
scale, we always quote our values for n s at a scale 
fco = 0.015/Mpc, where the tilt and running are 
less correlated, such that n s (ko = 0.015/Mpc) = 
n s (k = 0.002/Mpc) + ln(0.015/0.002)dn s /dlnA: [82[- 
An example of the remaining correlation between the 
spectral index and its running is shown in Fig. [3j 

III. RESULTS 

We now explore the constraints on extended pa- 
rameter spaces with the CMB (WMAP7+SPT), BAO 
distances (SDSS+2dFGRS), and an HST prior on the 
Hubble constant. Beginning with Sec. IIII CI we al- 
ways also consider SN distance measurements from 
the Union2 compilation. In Sec. IIII El we discuss the 
expected constraints from the Planck experiment. 



A. ACDM with Massive Neutrinos 

1. Enforcing the inflation prior on {ilk, dn s /d In k} 

In Table HU we begin by allowing the effective num- 
ber of neutrinos, sum of neutrino masses, and pri- 
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G(Q /0.25) L 

8 m 

FIG. 1. Joint two-dimensional marginalized con- 
straints on (T8(fim/0.25) ' 47 against {N e s, m„}. 
The black confidence regions (inner 68%, outer 
95%) are for the extended parameter combination 
"vanilla+A^gf-l- ^ m v +w+flk+dn s /d\n k" using the data 
from "WMAP7+SPT+# +BAO+SNe," while the ver- 
tical red lines denote the 68% confidence interval about 
the mean from the local (0.025 < z < 0.25) galaxy cluster 
abundance measurement of Vikhlinin et al. (2009) [891 ]. 



mordial helium abundance to vary as free parameters, 
both separately and jointly, in a ACDM universe. 

For ACDM alone, then with N e g and Y„ added sep- 
arately, we reproduce the results in Ref. In par- 
ticular, with TVeff = 3.87 ± 0.42 in the space given by 
"vanilla+A'eff ," we recover the reported 2a deviation 
from canonical N e g = 3.046 [13, El- Given the well 
known degeneracy between N e g and Y p [l(], HH E3] 
(also see discussion in SeclJ), we find N e g = 3.70±0.54 
when further allowing Y p to vary as a free parame- 
ter irrespectively of the BBN expectation. Allowing 
for three active neutrinos to have mass, and treat- 
ing additional contributions to N c g as masslcss, we 
find an even larger deviation with the standard value 
as N e g = 4.00 ± 0.43 for the parameter combination 
"vanilla-|-A^ e ff-|-^ m^" (consistent with [ll|). Here, 



the upper bound on the sum of neutrino masses is 
0.67 cV (95% CL) and we find the spectral index to 
be consistent with unity at la (n s = 0.987 ± 0.013). 
The neutrino mass constraint is to be compared with 
0.45 eV at 95% CL in "vanilla+J] m v " (consistent 
with [T3| and competitive with the robust upper 
bound at 0.36 eV when including CMASS 0). 

When we consider " vanilla+A/ r c fi'+X] m„+Y p " the 
statistical significance of the N e g deviation is reduced 
from 2.2a to 1.6a, and the upper bound on the sum 
of neutrino masses moderately weakens to 0.73 eV 
(95% CL). While the primordial helium abundance 
from the CMB+BAO+# has been found mildly 
in tension (~ 2a) [13, E| with that from observa- 
tions of metal-poor extragalactic H II regions [831 - 
HU, we find constraints on Y p consistent to within la 
with these observations. This is mainly due to the 
strong negative correlation between Y p and N c g (as 
reported in [13, [H El and detailed in Sec. HJ). For 
instance, Aver, Olive, & Skillman (2011) [83[ deter- 
mine Y p = 0.2534 ± 0.0083 via an MCMC analysis 
that accounts for both statistical and systematic un- 
certainties, which agrees with Y p = 0.277 ± 0.037 in 
"vanilla+Ar cff +r p " and with Y p = 0.261 ± 0.039 in 
"vanilla+A^eff+I] m„+Y p ." 

Moreover, when N e g and ^ m v are analyzed in a 
joint setting, we find that the data is both consistent 
with higher values of il c h 2 = 0.13 ± 0.01 and lower 
values of as = 0.80 ± 0.05, which perfectly agrees 
with low-rcdshift measurements of as from the abun- 
dance of clusters .891-1931] (as also noted in Ref. [HI). 
This is because the amount of suppression in mat- 
ter clustering by the free-streaming of light neutrinos 
increases with mass [U, [94T - I961 ] . which gives a large 
anti-correlation between ^2 m v and as , an example of 
which can be seen in Fig. [T] 



2. Relaxing the inflation prior on {Q.k, dn s /d In k} 

Let us now relax the strong inflation prior on the 
curvature of the universe and running of the spec- 
tral index by considering the parameter combination 
"v&ni\\&+N c fi+Y J ^h;+dn s /d\rik+VL k " in Table Iffll 

Here, N e g becomes increasingly consistent with the 
canonical value at 1.2cr (down from 2.2cr), mainly as 
a result of the anti-correlation with dn s /d\nk, which 
also brings the tilt down to n s — 0.978 ± 0.015 (from 
n s = 0.987 ±0.013). Further, we find that the correla- 
tion between ^2 m v and fifc degrades the upper bound 
on the sum of neutrino masses by close to a factor of 2 
to m v < 1-2 eV (95% CL). As a consequence of the 
well known anti-correlation with the sum of neutrino 
masses, which increases when relaxing the strong in- 
flation prior, the amplitude of matter fluctuations is 
seen to prefer smaller values at a 8 = 0.768 ± 0.070 (as 
compared to as = 0.798 ± 0.053 when dn s /dlnk and 
arc held fixed). 

However, as the strong inflation prior is relaxed, 
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■ WMAP+SPT+BAO+Jf +SNe 




FIG. 2. Joint two-dimensional marginalized constraints (inner 68% CL, outer 95% CL) for the extended param- 
eter combination "vanilla+iV fi-|- ^ m v -\-w+Q.k-\-dn s / d\n k," showing N e g against {w, Q c h 2 , n s } and 5^ m " against 
{w,fL c h 2 ,N eB }. The purple confidence regions are for " WMAP+SPT+ffo+BAO" and the black regions are for 
"WMAP+SPT+#o+BAO+SNe," where the BAOs and SNe are from SDSS+2dF and Union2, respectively. For Planck 
(T, E, (j>) in dashed white, we have centered the lcr error ellipses on the fiducial values of the parameters that went into 
computing the Fisher matrix. The exception to this convention is ^m, y , which we have shifted down to eV from its 
fiducial value of 0.17 eV for simpler visual comparison with the upper bounds from present data. 



both the running and curvature are consistent with 
zero to lcr. It is therefore far from certain that shifts 
in parameters other than {fife, dn s /dhik} that come 
about from relaxing the strong inflation prior are true 
manifestations that will hold with improved data. 



B. uiCDM with Massive Neutrinos 

1. Enforcing the inflation prior on {ilk, dn a /d\a k} 

In the previous section, we considered cases with 
neutrinos as massless and cases with neutrinos as mas- 
sive. However, as it is well established that neutrinos 
are indeed massive [28l - l3lj ]. we account for the sum 
of neutrino masses as a free parameter in all further 
treatments of neutrinos. In Table Hill we explore pos- 
sible degeneracies between {A^ff , ^ m^} and a con- 
stant EOS of the dark energy (w ^ —1). 

Beginning with " vanilla+w," we constrain a con- 
stant dark energy EOS: w = —1.10 ± 0.11 (as com- 
pared tow= -1.10 ± 0.14 without SPT). Consider- 
ing w in conjunction with "vanilla+TVofj-l- ^ to^," we 



find a reduction in N e s = 3.59 ± 0.57 (down from 
N e g = 4.00 ± 0.43), rendering it consistent with the 
canonical value to within lcr. This is caused by the 
w — N e ff correlation discussed in Sec Q] and shown in 
Fig. O Expectcdly, we also find a correlation between 
the dark energy EOS and the sum of neutrino masses 
(discussed in Sec. UJ, the latter of which degrades by 
close to a factor of 2 to X)m„ < 1.2 eV (95% CL). 

The joint impact of {N e g,J2 m v} on the dark en- 
ergy EOS is to weaken the constraint on it by roughly 
a factor of 3 to w = — 1.31 ± 0.30. Moreover, with 
the introduction of u>, the amplitude of linear matter 
fluctuations is mildly shifted to smaller values at erg = 
0.775±0.063 (compared to <j 8 = 0.798±0.053) because 
of the anti-correlation between w and erg that mainly 
enters through the growth function (e.g. see 0]). The 
spectral index shifts further away from unity to n s = 
0.968 ± 0.022 (down from n s = 0.987 ± 0.013). 

2. Relaxing the inflation prior on {Q.u, dn 3 /dlnk} 

In Table HLT1 we now consider the parameter combi- 
nation " vanilla+ N e ff+ ^2 77i„+ui" in conjunction with 



7 



■ WMAP+SPT+BAO+# +SNe 

WMAP+SPT+BAO+H 
□ Planck [T, E, if] 




FIG. 3. Same as Fig. [2] but for fife against {N e ff,^2m v ,dn s /d\nk} and dn s /d\nk against {iV c fi , m„, n s }. 



the running of the spectral index and curvature of the 
universe. We also consider a case with the primordial 
helium abundance as a free parameter. 

Given that we already identified separate degenera- 
cies between N e g — dn s /dhik and N c h — w in past sec- 
tions, it is not surprising that we obtain N e s = 3.10 ± 
0.74 to be in even closer agreement with the canoni- 
cal value for our extended parameter space. This is a 
result of the even more negative values preferred by 
dn s /d\nk = -0.018 ± 0.019 and w = -1.46 ± 0.39, 
shown in Figs. [5] and [3] However, the upper bound 
on J2 m f < 1-2 °V is robust to the further expan- 
sion of the parameter space, such that this bound 
holds for all three cases: "vanilla+A^ff-l- ^ m^+iu" , 
v vaml\a.+N e ff+J2 m u+^k+dn s /dlnk v , as well as 
" vanilla+TVeff-l- ^ m,y+w+^k+dn s /d In k" . 

In addition, when allowing for Y p as an in- 
dependent parameter (i.e. considering the case 
" vanilla+TVeff + J2 m u +w+Q, k +dn s / d In k+Y p " ), the 
upper bound on the sum of neutrino masses is only 
mildly weakened, while the error bars are large enough 
that the effective number of neutrinos is consistent 
with values of both 3 and 4 (to within 68% CL). In all 
of the above cases we continue to find erg consistent 
with that from cluster abundance measurements to 
within 68% CL, while the spectral index differs from 
unity by less than 95% CL. 

While our results are based on the construction of 3 
massive neutrinos, and (iV c ff — 3) massless degrees of 



freedom (as discussed in Sec. [TTJ) , we also considered 
imposing a hard prior of N c g > 3 in a new run with 
"va.ml\a,+N e ff+J2 m v+ w +Qk+dn s /dlnk", as this is 
the case where N e g would be the most affected by the 
prior. Given the weak correlation between ^ m„ and 
iVeff, the upper bound on the sum of neutrino masses 
doesn't change, while the data is still consistent with 
no extra relativistic species as N Q s = 3.65 ^qq' ^qq, 
where the two sets of upper and lower boundaries de- 
note 68% and 95% CLs, respectively. It is clear that 
our findings are qualitatively unchanged with this al- 
ternative choice of prior on the effective number of 
neutrinos (also see captions of Tables IIIIHIV1 and lVip . 



C. miCDM with Massive Neutrinos, Running, 
and Curvature: Including Supernovae 

Since much of the work in bringing N c s in agree- 
ment with the canonical value is done by the possibil- 
ity of evolving dark energy, for which the constraints 
from the CMB, Ho, and BAO measurements that we 
have considered are relatively weak, we further include 
SN data from the Union2 compilation in order to more 
effectively constrain the dark energy EOS and param- 
eters with which it strongly correlates. 

In Table [TV] we find that the addition of SN ob- 
servations help constrain the dark energy EOS to 
w = — 1.05±0.07 when analyzed along with the vanilla 
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TABLE V. Constraints on Cosmological Parameters using SPT±WMAP±# ±BAO±SNe. 







io(a)CDM 


m(a)CDM 


w(a)CDM 

rf. In 


(u(a)CDM 

1 n A: 


_T 1 111 lOyl V 


100fl b h 2 


2.226 ±0.042 


9 94Q -1- O 0/17 








lOOQc/i 2 


11.37 ±0.46 


13.3 ±1.1 


13.3 ± 1.1 


14.1 ± 1.4 




10 4 6» 3 


104.11 ±0.15 


103 96 + 18 

lUv, J VJ _l_ U. -LU 


103 99 + 18 

luU. Jt7 _1_ VJ • J- (J 


103 66 + 36 




T 


0.083 ±0.014 


VJ . UOU _l_ U.uU 


n (189 + 01 6 


n nqn + 01 6 

VJ . UC/L/ _1_ VJ . VJ i VJ 




Its 


0.963 ±0.011 


978 ±0015 

Wtl/ 1 U _1_ VJ . VJ J- <J 


970 + 019 


950 + 026 

VJ ■ <J \J\J 1 VJ . \J£i VJ 




In (10 10 A a ) 3.203 ±0.041 


3 1 94 -\- 044 


3 1 94 ± 054 


q 200 + 055 

O • VJ U _l_ VJ . VJ i J !J 


Extended 


w 


-1.10 ±0.17 


-1.05 ± 0.19 


-1.08 ±0.21 


-1.12 ± 0.21 




w a 


0.20 ±0.64 


-0.4 ± 1.0 


-0.3 ± 1.2 


-0.3 ± 1.2 




N cS 




3.84 ±0.45 


3.57 ±0.59 


3.75 ±0.68 




E» [eV] 




< 1.2 


< 1.4 


< 1.8 




d In 






-0.012 ±0.020 


-0.038 ± 0.030 




ioon fc 






0.7 ±1.1 


1.3 ±1.2 












0.168 ±0.079 


Derived 


^8 


0.827 ±0.046 


0.779 ±0.062 


0.763 ±0.076 


0.736 ± 0.084 



Same as Table ITVl but for a time-dependent parameterization of the dark energy equation of state, of the form w(a) = uiq + 
(1 — a)w a (as opposed to time-independent w). 



parameters (35% reduction in uncertainty compared 
to no SNe). This constraint degrades to —1.10 ± 0.11 
when expanding the parameter space to further in- 
clude {N c g , m v , f2fc, dn s /dhi fc}, but is still a fac- 
tor of 4 stronger than the equivalent case where SNe 
are not included in the analysis. Improving the con- 
straint on w is helpful in breaking much of the de- 
generacy between dark energy and the effective num- 
ber of neutrinos, resulting in N e g = 3.88 ± 0.44 for 
the case of " vanilla± N e g+J2 m u +w" (as compared to 
N e g = 3.59 ± 0.57 without SNe, and as compared to 
N e g = 4.00 ± 0.43 with a prior w = —1). However, 
as before, when relaxing the strong inflation prior on 
{flk, dn s /d In k} we find the effective number of neu- 
trinos becomes consistent with the canonical value to 
68% CL (as N cS = 3.58 ± 0.60). 

Low-redshift SN measurements are useful in re- 
ducing the correlation between {J^ m„, w, Ho}, which 
drives the 1.2 eV (95% CL) upper bound on the 
sum on neutrino masses for the case "vanilla±./V e ff + 
Y^rriv+w" down to 0.9 eV (Tables HQ and |HT|) . 
However, since SN observations do not much im- 
prove the constraint on the curvature when added 
to CMB+#o+BAO (shown in Fig. EJ), the SNe are 
unable to lower the upper bound on the sum of 
neutrino masses from 1.2 eV (95% CL) when re- 
laxing the strong inflation prior (i.e. for the case 
"vanilla±./V e ff±^TO I/ ±w±17fc±dn s /dlnfc"). The pa- 
rameter that most strongly increases the upper bound 
on the sum of neutrino masses when singularly added 
to "vanilla+^TOy" is the curvature, which renders 
X>„ < 1.0 eV (95% CL) in "vanilla+^m^+Ofc." 

Expanding the parameter space to allow Y p to 
vary as an independent parameter (i.e. consider- 
ing "vanilla+A r ff± ^m^+w+Qk+dns/dlnk+Yp"), 



we find a mild shift in N e g = 3.78 ± 0.61 (as com- 
pared to N c s = 3.58 ± 0.60), and a stronger shift in 
Y^, m v < 1-7 eV (as compared to X) TO f < 1-2 eV at 
95% CL). Meanwhile, Y p = 0.176±0.079 shows a pref- 
erence for lower values but is still consistent with mea- 
surements of Y p from low-metallicity H II regions (83 - 
|88j. For all of the non- minimal cases considered in 
Table IIV1 n s is consistent with unity to at least 95% 
CL, and as mildly prefers values less than 0.8 but is 
still greatly consistent with cluster abundance mea- 
surements as shown in Fig. [TJ Moreover, we find that 
larger values of the dark matter density are preferred, 
as f2 c /i 2 generally lives around 0.13 ± 0.01. 

For the particular parameter combination that 
shifts N e g the closest to a value of 3 from above 
(i.e. "vanilla±7V e ff±^m ly +i«+f7fc±dn s /(ilnfc"), we 
also considered a run with the prior N c g > 3 imposed. 
Here, we continue to find the effective number of neu- 
trino species to be consistent with the standard value, 
as JVeff = 3.74 3 qq : 3 oo, where the two sets of upper 
and lower boundaries denote 68% and 95% CLs, re- 
spectively. The constraints on other parameters such 
as w and ^2 m v change by less than 10% with this 
alternative choice of prior. 

Next, we move on to other parameterizations of the 
dark energy, such as the popular expansion w(a) = 
wq + (1 — a)w a and an early dark energy model in 
which the EOS of the dark energy tracks the EOS of 
the dominant component in the universe. 



D. Alternative Dark Energy Parameterizations 

Given our ignorance of the nature of dark energy, 
once we move away from a cosmological constant, 
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TABLE VI. Constraints on Cosmological Parameters using SPT+WMAP+H +BAO+SNe. 







eCDM 


eCDM 


eCDM 


eCDM 
+N cft +Y,m IJ +Y p 

± J W.\ + ^fc 
In k. 


L I llllell V 


1 DflOi./) 2 


ooq _)_ n 941 


9 nr\c\ 4. n 047 




9 1 70 _j_ n nsi 




1000 h 2 


ii Kf: _i_ n 40 

ll.JU _1_ L/.rlii 




13.4 -J-_ 1.1 


13.9 ±1.5 




10 4 6* s 


in4 no + o 17 

-L\JT: a \J\J _1_ U. 1 1 


103 89 + 18 


iqq on _|_ n 1 q 


1 03 72 + 39 




- 


0.085 ± 0.014 


0.090 ± 0.015 


0.091 ± 0.015 


0.094 ± 0.016 






966 + 01 1 


n 984 + 015 


966 + 019 


957 + 027 




in nn 1( M i 




q 1 74 4. n 942 


q 1 7fj 4- n nc;i 

u>l 1 U _1_ U.UJ1 


q 1 «2 + 055 


Extended 




-1.082 ±0.079 


-1.11 ±0.10 


-1.16±0.12 


-1.17 ±0.14 




n, : 


< 0.030 


< 0.025 


< 0.049 


< 0.049 




N efl 




3.85 ± 0.43 


3.24 ±0.63 


3.37 ±0.68 




■£m„ [eV] 




< 0.96 


< 1.6 


< 2.0 




d In k 






-0.023 ±0.021 


-0.034 ±0.032 




won k 






1.5 ± 1.2 


1.7 ± 1.3 












0.206 ± 0.086 


Derived 


C8 


0.805 ± 0.045 


0.773 ± 0.063 


0.703 ±0.095 


0.692 ± 0.097 



Same as Table ITVl but for an early dark energy model with present EOS too an d density at high redshift Q e (as opposed to time- 
independent to). We report the 95% upper limit on Q e (and ra v as before). For the "eCDM+iVeff + E m v + dn s /d In k ± Q^" 
case, we also considered a run where we impose a hard prior of iV e fj > 3. Here, we find N e g = 3.60 3 qq' 3 'oo' where the two 
sets of upper and lower boundaries denote 68% and 95% CLs, respectively. The largest changes this prior induces in other 
parameters are seen in n a = 0.975 ± 0.015 (compared to n s = 0.966 ± 0.019), dn s /d\nk = —0.014 ± 0.017 (compared to 
dn a /d\nk = -0.023 ± 0.021), lOOCfe = 1.0 ±1.1 (compared to 100f2 fe = 1.5 ±1.2), and U e < 0.042 (as compared to C e < 0.049). 
All of the other parameters are modestly affected (< 10%) by our choice of prior on N c ff. For the vanilla±too±^e case, we al so 
considered a run with a hard prior w > —1, for which we find Q e < 0.023 at 95% CL (as compared to Q e < 0.019 in Ref. |l06l |). 



there is no adequate reason to restrict our analyses 
to a constant EOS from the point of view of parti- 
cle physics, in particular if we wish to describe the 
dark energy as a scalar field or modification of grav- 
ity (e.g. (H, |H, H3-H9|). Thus, as an extension of 
the previous section, we now consider models of the 
dark energy in which the EOS varies with time. While 
SN measurements proved useful in breaking parameter 
degeneracies with a constant EOS, we aim to under- 
stand how well these degeneracies are broken for less 
constrained dark energy models. 



1. Late-Time Dark Energy with Evolving 
Equation of State 

The first of our alternative parameterizations for 
the dark ener gy is given by the two-parameter 
model [99MT0ll] adv ocat ed in the report of the Dark 
Energy Task Force flO^ : 



w(a) = w + (1 - a)w a , 



(3) 



where wq is the EOS at present, while conventionally 
w a = — 2dw/dlna\ a=1 /2 |99l llOOt ] . Eqn.[3]may also be 
viewed as a first order Taylor expansion of the EOS, 
where w a = —dw/da\ a =i- 

As compared to the case with a constant dark en- 
ergy EOS, we find that the new parameterization for 
late-time dark energy doesn't significantly change our 



constraints on other cosmological parameters. Ex- 
pectedly, the two most sensitive parameters are fi^ 
and ^to„, the constraints on which degrade by less 
than 15% and 30%, respectively. 

In Table fVl we constrain wq = —1.10 ± 0.17 and 
w a = 0.20 ±0.64 for the extension "vanilla+tuo+uv." 
When we in stead use SNe from the "Constitution" 
compilation |l03j |. we find wo = —0.93 ± 0.13 and 
w a — —0.36 ± 0.65, which are consistent with the 
constraints on these parameters in Ref. Q (perfect 
agreement when excluding SPT). Clearly, precise SN 
measurements are critical to understanding the true 
values of these EOS parameters. In further extensions 
of parameter space, the constraint on wo degrades by 
up to 20%, while the constraint on w a degrades by up 
to a factor of 2. 



2. Early Dark Energy 

Late-time dark energy models suffer from the well 
known coincidence problem. The value of the dark en- 
ergy density has to be fine-tuned so that it only affects 
the dynamics of the universe at present. This coinci- 
dence problem motivates the exploration of models in 
which the evolution of the dark energy density is such 
that it is large enough to affect the universal dynamics 
even at z > 2. 

A realization of early dark energy (EDE) is given 
by the "tracker" parameterization of Doran & Rob- 
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WMAP+SPT+BAO+H +SNe 
Planck [T, E,c|>] 




FIG. 4. Joint two-dimensional marginalized constraints on the early dark energy density Sl e against {N e a, Y rn u } for the 
extended parameter combination "vanilla+iV e fj+ Yl m v -\-wo-\-Q. e -\-Q,k+dn s /d\D. fc." The black confidence regions (inner 
68%, outer 95%) are for " WMAP+SPT+# +BAO+SNe," while the forecasted la error ellipses for Planck temperature, 
E-mode polarization, and lensing potential power spectra (T, E, (/>) are shown in dashed red. Although the Fisher matrix 
constraints on the parameters {^m„,f2 e } were evaluated at {0.17 eV,0.01}, they have been shifted down to {0,0} for 
simpler visual comparison with the upper bounds from present data. 



bers (2006) [l04| . where the dark energy tracks the 
dominant component in the universe. In a sense, it is 
simpler to parameterize the dark energy density evo- 
lution directly, rather than express it in terms of an 
evolving equation of state. We use a modified form of 
the original parameterization th at tracks t he equation 
of state of the dominant energy 

EMEU, 

n d (z)-n d0 hl{z) 

( (14- z ) 3+3w ° \ 
+ (l-L_J_), (4 ) 

h 2 w (z) = n M (i + z) 3+3wo + o ro (i + zf 
+ a r (i + z) 4 + n k (i + z ) 2 , 

where {O r , £l m } are the present radiation and matter 
densities in units of the critical density. The present 
matter density is further composed of the densities of 
the cold dark matter, baryons, and massive neutrinos 
(O m = f2 c + Qf, + Q v ). 

As described in Ref. [l05j . the function v(z) should 
have the properties that it asymptotes to unity at 
large rcdshift and u(0) = 0, thus ensuring that fld(z) 
asymptotes to f2 e at large re dshift and f2d(0) = &do- 
We use v(z) = 1 - (1 + z) 3w ° [Toi|, but any other pa- 
rameterization such that d\u(v) / dln(z) = 0(1) will 
give similar results. Note that the first term propor- 
tional to firfo is the dark energy density as a function 
of redshift for a model with present density of dark en- 
ergy Qdo and constant EOS wq. Thus, in this param- 
eterization with early dark energy, the effect at low 
redshift is the same as a model with constant EOS. 

By approximating the effect of a dark energy com- 
ponent with time- varying EOS using the PPF module 
of Ref. [z3| , we allow wq to freely vary above and be- 
low the w = — 1 boundary, unlike the treatments in 
Rcfs. [2~il . 1 106j | . We compute the equation of state us- 



ing the expression w(z) = — 1 + din ^ d ^ 1 ™^ . 
At z=0, w(z) = wq and increases with z, tending to 
when the dominant component of energy density is 
due to pressureless matter, and to 1/3 when the uni- 
verse is dominated by radiation. Quantitatively, the 
O e term ( "early dark energy" ) in Eqn. 2] constitutes 
[0,2.1,8.0,17.7]%, at redshifts z = [0,1,2,3] respec- 
tively, of the overall amount of dark energy Sld(z) for 
w = -l and fl e = 0.01. 

The impact of early dark energy on the consid- 
ered observables mainly comes through increasing the 
expansion rate and in shiftin g the ma tter-radiation 
equality to a later epoch [24], Il04| - |l06| . We can use 
the CMB to constrain an EDE model via its effects 
on {z eq , S , 9d}, while the improvement from BAO and 
SN distance measurements are modest, as the expan- 
sion rate in a model with EDE is designed to mas- 
querad e that of models with late-time (z < 1) dark 
energy |l04lll05j . For our EDE model we fix the sound 
speed c s = 1 , whi le other choices have been explored 
in Rcfs. pMH07t. 

For the " vanilla-|-iuo+^e" case, we find f2 e < 0.030 
(95% CL). We also considered a run with a hard prior 
w > -1, for which we find Q e < 0.023 at 95% CL 
(as compared to f2 e < 0.019 in Ref. |l06| ). When ex- 
panding the parameter space to include the neutrino 
sector, we find a 20% reduction in the upper bound 
on Q e , due to its correlation with another one-tailed 
distribution ^2 m v ■ As compared to the case where 
the EOS is described by a simple constant, we find 
modest changes in the constraints on all other param- 
eters. However, as a result of the correlations between 
f2 e and {fifc, dn s /d\n k}, we find N c g = 3.24 ± 0.63, 

m u < 1-6 cV, and erg = 0.703 ± 0.095 when relax- 
ing the strong inflation prior (as compared to N e g = 

3.58 ± 0.60, m » < !- 2 eV > and °"8 = °- 774 ± °- 072 
when w is a constant). The upper bound on the early 
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dark energy density itself degrades by a factor of 2 
to Vl e < 0.049 (95% CL). To obtain these shifts, the 
curvature and running show weak (lcr) preferences for 
nonzero values. 

For the same parameter space, we also carried 
out a run with an explicit N e s > 3 prior, finding 



CMB Survey Properties 



Experiment Channel FWHM AT/T X 10 6 AP/T X 10 6 



N, 



elf 



3.60 



3.74, 4.47 
3.00, 3.00 



, where the two sets of upper 



and lower boundaries denote 68% and 95% CLs, re- 
spectively. This prior lowers the upper bound on the 
EDE density to fi e < 0.042 at 95% CL (as compared 
to f2 e < 0.049), while the constraints on other pa- 
rameters such as w and YJ m v change by less than 
10%. When further including y p as a free parame- 
ter, we find qualitatively modest changes in our con- 
straints, similar in nature to those discussed in sec- 
tions [ml] and unci 



E. Parameter Forecasts for Planck 

Having discussed the present status of constraints 
on expanded parameter spaces with CMB, Hq, BAO, 
and SN measurements, we next explore the constraints 
from CMB temperature, E-modc polarization, and 
lensing potential power spectrum measurements with 
Planck jl08l . Il09| . To this e nd, we employed a Fisher 
matrix formalism jl ldL llll| , such that the parameter 
covariance matrix is given by the inverse of 



^ — \j n 



dp a dpp 



(5) 



where the CMB temperature (T), E-mode polariza- 
tion (E), and lensing potential (4>) power spectra enter 
the symmetric matrix 



/ C M C 4>T Q 

c T <t> r<TT r<TE 

W °£ W 

V o Cf T Cf E 



(0) 



power spectra contribute additively to 
/-y 2 (2/(2^+l)) 1/2 {C ah (t) + 5 ab N ab {£)), 



The noise 

c ab (£) = , sky 

where {a, 6} 6 {</>, T, E}. The experimental specifica- 
tions are listed in Table I VIII For the terms in Eqn. [SJ 
we carried out two-sided numerical derivatives with 
steps of 2% in most parameter values. We have con- 
firmed the robustness of our results to other choices 
of step size. For further details on our prescription, 
inclu ding how to obtain the noise power spectra, see 
Ref. jl05| . Therein, we also consider weak lensing 
tomography, galaxy tomography, supernovae, and ex- 
tensive set of cross-correlations for future wide and 
deep surveys. 

As shown in Figs. [21 [31 [H Planck will be extremely 
helpful in improving the constraints on extended pa- 
rameter spaces. At the la level, considering the com- 
bination "vanilla+A' c ff+ YJ m l ^+w+Qk+dn s /dhik , \ 
the effective number of neutrinos could be constrained 



Planck 


100 


10 


25 


10 




143 


7.1 


16 


30 




217 


5.0 


21 


19 



TABLE VII. Experimental specifications for the Planck 
mission. The sky fraction / s k y = 0.65, and the angular 
multipoles extend from £ m i n = 2 to ^ max = 2000. The 
channel frequencies are given in GHz, and the angular res- 
olutions in arcminutes. 



to a(N cff ) = 0.23, mainly from the CMB tem- 
perature power spectrum, and the sum of neutrino 
masses to a(J2 m v) = 0-19 e V, mainly from the 
CMB lensi ng potential power spectrum (consistent 
with Refs. [To! tm El)- These constraints on N cS 
and YJ m v are a factor of 3 stronger for both parame- 
ters than the present constraints from a joint analysis 
of "WMAP+SPT+iJo+BAO+SN" data. 

When further allowing for the possible exis- 
tence of a non-negligible component of dark en- 
ergy in the high-redshift universe (i.e. consider- 
ing "vanilla-l-A'eff-r- YJ m^+WQ+fl e +flk+dn s /d In fc"), 
Planck constraints on VJm„ and N e g only degrade 
by 10% and 20%, respectively. This is because the 
CMB temperature power spectrum achieves a factor 
of 6 improvement in the constraint on the EDE den- 
sity at a(Q e ) = 0.0087 (which only improves by 5% 
in the full analysis) , removing much of the degeneracy 
with other parameters obtained from the CMB. Ex- 
pectedly, the orientation of error ellipses for Planck 
and present data in Figs. [21 El HI match for parame- 
ters that are mainly constrained by the CMB (T,E), 
while they differ for parameters, such as w, for which 
the CMB temperature and E-mode polarization pro- 
vide inferior constraints. 

Moreover, we allow for the primordial helium 
abundance to vary in " vanilla+A^fj-l- YJ m^+Vp" , 
where the la constraints on {N D g , YJ to„, Y p } are at 
the level of {0.25, 0.14eV, 0.015}. In the extended 
space "vanilla-r-Af e ff-|- YJ m„+w+flk+dn s / din k+Y p ," 
these constraints degrade to {0.27, 0.20eV, 0.019}, 
respectively (such that these neutrino constraints 
are comparable with the Q e case). However, the 
expected parameter constraints worsen signifi- 
cantly when considering {N c h , YJ m„, Y p , fi e } in 
conjunction. For the maximal parameter extension 
('Wa,mll&+N eS +J2m v +Wo+n e +Qk+dn s /dlrik+Y p "), 
we find la constraints on these four correlated pa- 
rameters at the level of {0.51, 0.29eV, 0.030, 0.013}, 
respectively. In this scenario, Planck will need to 
be combined with external datasets to break the 
parameter degeneracies, in particular low-redshift 
measurements of the expansion history. 

Thus, assuming a strong BBN prior on Y p (alter- 
natively, at less than 10% degradation, assuming the 
dark energy is a pure late-time phenomenon at the 
level of fi e < 3 x 10~ 3 ), it is expected that Planck 
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alone will be able to determine the possible existence 
of extra relativistic species to 4cr confidence and the 
sum of neutrino masses to 0.2 eV, regardless of the 
extent of the remaining parameter space. With the 
ability to strongly constrain ^m„, there is promise 
for Planck to find evidence for nonzero neutrino mass, 
in particular w hen c ombined with prob es of the large- 
scale structure [l05] (also see [ll4fl23j ). 



IV. CONCLUSIONS 

With the latest cosmological data sets of the 
cosmic microwave background (WMAP7+SPT), 
baryon acoustic oscillations (SDSS+2dF), supernovae 
(Union2), and the Hubble constant (HST), we have 
explored in closer detail the dependence of constraints 
on the effective number of neutrino species and the 
sum of neutrino masses on our assumptions of other 
cosmological parameters, including the curvature of 
the universe, running of the spectral index, primor- 
dial helium abundance, evolving late-time dark en- 
ergy, and early dark energy. 

In a combined analysis of the effective number of 
neutrinos and sum of neutrino masses (with 6 other 
ACDM parameters), we find mild (2. 2c) preference 
for additional light degrees of freedom. However, the 
effective number of neutrinos is consistent with three 
massive neutrinos and no extra relativistic species to 
la when further including {w, f2^, dn s /dlnk} as free 
parameters. The transformation of a constant EOS 
to one that varies with time (wo,w a ) doesn't signifi- 
cantly change the constraints on N e g and ^2 m v (less 
than 10% and 30%, respectively). The agreement with 
JV e fF = 3.046 improves with the possibility of an early 
dark energy component, itself constrained to be less 
than 5% of the critical density (95% CL) in our ex- 
panded parameter space. 

Added to a minimal ACDM universe, m v < 
0.45 eV (95% CL), dominantly from WMAP+HST. 
The sum of neutrino masses is bounded at 1.2 eV (95% 
CL) when jointly allowing for a constant dark energy 



equation of state, curvature, and running to vary as 
free parameters. The upper bound degrades to 2.0 eV 
(95% CL) when further including the primordial he- 
lium abundance and early dark energy density as addi- 
tional degrees of freedom. The single parameter that 
most strongly increases the upper bound on the sum 
of neutrino masses when added to "vanilla+^ rn„" 
is the curvature of the universe, which weakens the 
bound by more than a factor of 2 to m v < 1-0 eV 
(95% CL) in "vanilla+X^+fifc." 

In extensions of the standard cosmological model, 
the derived amplitude of linear matter fluctuations 
as is found consistent with low-redshift cluster abun- 
dance measurements to within la, and the spectral 
index agrees with unity to within 1 to 2 a. Moreover, 
larger values of the dark matter density are preferred, 
as il c h 2 generally lives around 0.13 ± 0.01. When al- 
lowing the primordial helium abundance to vary as 
a free parameter, we consistently find la agreements 
with estimates from observations of low-metallicity ex- 
tragalactic H II regions. 

With the advent of increasi ngly sensitive CMB and 
large-scale structure data [HI, Il08l our 
understanding of the neutrino sector depends criti- 
cally on the ability to distinguish its signatures from 
other cosmological parameters. Fortunately, even for 
extended parameter spaces, Planck alone could deter- 
mine the possible existence of extra relativistic species 
at the 4<T level and constrain the sum of neutrino 
masses to 0.2 eV (68% CL). Next-generation probes 
of the low-redshift expansion history and large-scale 
structure hold the key to further improving these es- 
timates. 
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